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ABSTRACT 


A study of an Automatic Power Compensator with inputs of airspeed 
and angle of attack was conducted by means of a simplified analytical 
analysis and analogue simulation. Digital Computer programming was 
used for data reduction. The Chance Vought F-8 aircraft was used as 
the vehicle for the study. Loop gains were initially set by means of a 
simplified analysis based on the eee phugoid and short period 
motions. The validity of such an analysis was demonstrated by ana- 
logue simulation of the complete system. The system is feasible: both 


for independent use or as a component of a complete automatic flight 


system. 
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TABLE OF SYMBOLS 

(For those symbols not defined in the text) 
Aerodynamic center 
(I) Moment of Inertia about the Y axis 
chord (MAC of wing) 
coefficient of Lift 
coefficient of Drag 
coefficient of Moment 
center of gravity 
acceleration of gravity 
angle of attack proportional gain 
airspeed proportional gain 
airspeed integral gain 
Lift force 
tail length - distance from cg to ac of horizontal tail 
mass 
vertical acceleration perturbation (g units) 
horizontal velocity perturbation 
steady state velocity 
reference velocity (V = U,) 
vertical velocity perturbation 
aircraft weight 


reference area (wing area) 


ll 


S Laplace operator 


T thrust 
(t) indicates variable in time domain 
oC angle of attack, a (see Fig 2) 


O 
elevator deflection (see Fig 2) 


3 


angle above horizontal (see Fig 2) 


@® 


flight path angle (see Fig 2) 


atmospheric density 


=) a) 


time constant for delays (see Table IV) 


indicates determinant expansion 


indicates matrix formulation 


Superscripts: 

d( ) 

dt 

marie 2) 

- dt2 
Subscripts: 
O initial or steady state value 
t error quantity or signal 
cS command signal 
f delayed signal 


Aerodynamic derivatives 


The standard NACA notation is used for the non-dimensional derivatives 


in that 
C C 
Clu = So, Cn~ = oe > etc. (See Table II) 
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L.. Introduction. 

This study is a continuation of a series of investigations designed 
to evaluate various types of "Auto-throttle” devices. The study was in- 
stituted in the Aeronautical Engineering Department of the United States 
Naval Postgraduate School in 1963. §, 9] 

The need for such a device on high performance carrier aircraft and 
the basic concepts relative to the system have been previously discussed. 
Briefly, the power compensator is a device meant to automatically stabi- 
lize the airspeed at a pre-selected value without attention from the pilot. 
The pilot then maintains the desired attitude or glide slope. This concept 
is not limited to, but is especially desirable during carrier type approaches 
where modern, high performance, swept-wing aircraft must operate ina 
range of inherent speed instability. With the heavier aircraft now being 
programmed for carrier use lower landing speeds are most desirable. By 
allowing the pilot to focus his attention on only one correction during ap- 
proach, .a-power compensator permits the safe utilization of slower approach 
speeds. Basically then, the automatic power compensation is a means of 
artificially tailoring the aircrafts inherent flight characteristics in the low 
speed region to meet the opposing requirements of safety and lower landing 
speeds. In this sense such a device could prove applicable; even neces- 
sary, to the landing of the larger super-sonic transport aircraft now under 
development. 

The initial study demonstrated the need for such a device by showing 


real time, pilot contrclled, analogue simulated mirror approaches with and 
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without an auto-throttle. [8 | The second paper attempted to improve the 
systems by using different inputs. [9] The systems analyzed thus far use 
inputs of airspeed/pitch angle and angle of attack/normal acceleration. 
These are systems of the type investigated by the United States Navy for 
operational use. The angle of attack/normal acceleration system is pre- 
sently in use on operational aircraft. 

Utilizing the results of the previous studies, Professor E. J. Andrews 
of the U. S. Naval Postgraduate School proposed a much simplified analytic 
means of setting loop gains. [4] 

This paper will investigate a power compensator that uses inputs of 
airspeed and angle of attack. Within the knowledge of the author no such 
system is in use at present. This type of a system might be thought of as 
a "hybrid system" since it utilizes one of the inputs from each of the two 
systems investigated thus far. The motivation for such a system lies in the 
fact that airspeed and angle of attack are two of the most "basic" governing 
parameters of flight. Since these parameters are related by 

(a+ mW = + P V2SCy 6% 
the values of airspeed and angle of attack will be related to the aircraft 
weight and the desired type of flight path. Fora given aircraft weight, the 
longitudinal flight path is determined completely by the elevator position and 
the throttle setting (airspeed). 

Therefore, with such an auto-throttle stablizing airspeed, the pilot 


(or autopilot) need only control the pitch attitude to obtain any desired flight 
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path. This type of an auto-throttle would be especially useful for flying 
a pre-determined flight path designed to give maximum range, maximum 
rate of climb or a desired glide path.over a range of airspeeds. 

In the investigation of this type of a power compensator, the simpli- 
fied analytic approach will be demonstrated in approximating the desired 
gain settings. The Chance-Vought F-8 aircraft was.used as the vehicle in 
the previous studies and will be used in this study. 

The equations of motion are solved using an analogue computer. A 
digital computer is used to determine aircraft parameters and roots of char- 
acteristic equations. 

All aerodynamic parameters are used in the English dimensional form. 
All equations and programming are given in real time. 

The aerodynamic data used was obtained from the Chance Vought 
Division of Ling-Temco Vought and an earlier report on an evaluation of a 
power compensator. [6, 7| 

A word of appreciation is due Professor E. J. Andrews for his assis- 
tance during this investigation. 

Ce Discussion. 
2.1 General 

The system under study is actually equivalent to the angle of attack/ 
normal acceleration system evaluated previously. The auto-throttle equa- 
tion for the </n system is: 


Eq(2-1) To = Kix - K 


gn 


thes 


1S 


Since L = 4 p(v'+a)28Cz, = ie SCL (Veith 2uV.+ UZ) 





alias MG? » OL ~ 1psc, (2v + 2u) £ SCV 
Then _ ire ts eal sa p(2V + 2u) LYP 
nW = AL 
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og 


On _ 70m, fo, . fe 
——— : alien W 
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OY 
c 
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Therefore: 


n= + PV2-B-Cr oc + p2-CrVvu = K30° - Kyu 


The Auto-throttle equation becomes: 


To = Koc - Kj (K3X + Kyu) = (Ky - K,K,)¢- K Kyu 


c 24 


Eq(2-2) To = Ke® - Kyu 


A block diagram of the complete system is shown in Figure 1. 


2.2 Airframe Equations 

The 3-degree of freedom, small perturbation, -aircraft longitudinal 
equations were used to solve for the airframe motion. The equations shown 
below are as developed in any standard text on aerodynamic stability. [3 ] 
The definitions of the derivatives are listed in Table I. The dimensional 
derivatives were computed from the non-dimensional derivatives utilizing 
a digital computer program given in the Appendix. The equations used for 
this conversion are given in Table II. The values of the dimensional deri- 
vatives obtained are listed in Table III along with all other pertinent flight 
parameters. 
Eqs (2-3) i = Xyu + Xyw + XgW - gcos@@ + XT + Xy/) 


(1 - Zy)W = Zu + Zyw + (Zy + Us) - gsin@ @ + ZpT + Zyl) 
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@ = Myu+ Myw+ Myw+ M0 + MgT + My7) 
The angular and velocity relationships are defined in Figure 2. The 
derivatives concerning thrust are dervied in Figure 3. These equations were 


reduced by inserting the values. 





and rearranged to give 


Eqs (2-4) u = Aju + U.X% - gcoos®,8 + XpT + X>!) 
= Thy + 2,04 § - -Tsind 9+ “Tp y Su 
, o Ww Uo Up Uo 


@ = Mut UpMyw% + UpMy,% + Mg@+ MyT + My’) 


2.3 Power Compensator 


The equation of the power compensator is: 


Eqs (2-5) To = pea T 
x 


K 


k 
Ky 
i ge ee 


where J~ and 7s are estimates of the angle of attack and airspeed measure- 
ment delays respectively. It can be seen from this equation that the power 
compensator will apply a forward (+) thrust for a nose up (+) angle of attack 
perturbation and an increase in airspeed (+ u,) will cause a reduction of 
thrust. The compensator utilizes angle of attack and airspeed proportional 
feedbacks plus an airspeed integration to insure zero steady state errors in 


the airspeed loop. 


ie 


2.4 Transfer Functions 
The transfer functions used are defined in Table IV. The delays were 
obtained from various sources. 6, 7| The airframe transfer functions are 
derived in the Appendix from the airframe Laplace equations shown below 
in matrix form. The airframe transfer functions were evaluated by means 


of a digital computer program given in the Appendix. 


Eqs (2-6) ay AG = Or + gcos®, u 
— + (S - Z,) _ gi —S8iueo ol 
U 
oO O 
UL 
“Ma  - Ug (My + SMz) sz-ms | | 6 | 
mr Xy 
T 
= ZT Zy 
Uo we 


7) 


Mr My | 
The open and closed loop transfer functions were derived from 


Figure 1 using standard methods. 





Open Loop: at ay al! 
Eqs (2-7) | 
“U = Di Ic u 
_ - i. 
Ue Ua T 
Eqs (2-8) 
ee : 
| Tee 'E u K 
ae aye oe 
Ug: Te Te T ce To TE T 


Only the airspeed transfer functions are given above. Transfer func- 


tions for ~/~Ke, u T) : aa could also be developed easily from the above 
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expression since the denominator, ie. the characteristic equation, would 


be the same in all cases. 


2.9 Digital Computer Programs 

The Digital Computer programs used in the previous studies were 
rewritten in FORTRAN 63 in a subroutine format. These programs are used 
to convert either wind tunnel data or non-dimensional derivatives to the 
dimensional derivatives necessary for the Analogue Program. These deriv- 
atives are then used as inputs to other subprograms which analize the basic 
airframe motion and compute the airframe transfer functions for the inputs 
of elevator and thrust. Thus it is possible to proceed directly from wind 
tunnel data to the response of the airframe for several flight conditions, 
using these results to analize any gain changes necessary in the control 


loops. These programs are contained in the Appendix. 


2.6 Complete Analysis 

The above transfer functions, both open and closed loop, could now 
be used to analyze the system. This involves a 5th order characteristic 
equation with at least three variable gains. The numerators are of third 
order orless. The characteristic equation is evaluated in the Appendix 
for later use in a Root Locus Diagram. At this point the methods of Root 
Locus plots combined with frequency response methods could be used to 
analyze the response for various combinations of gains to elevator inputs, 


and step gusts of uand &. This procedure would require an attack by 
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digital computor methods for at least computing the roots of the various 
polynomials generated for each value of the gains used and then perhaps 
estimating the response obtained by means of simplified Bode plots or 
Root Locus plots. But this rigorous, automatic control engineering appro- 
ach at this point is unnecessary and uneconomical if not also mathematic- 
ally unsound. Within the accuracy of the assumption inherent in the basic 
airframe equations (Eqs. 2-3) a much simplified approach is justified. In 
fact, it is known that the responses calculated by this "exact" method are 
not only in error, but for the case of shear gusts, can give unconservative 


results for structural loading parameters. [3] 


2.7 Simplified Analysis 
This analysis is based on the standard assumption of reducing the 
airframe equations of motion by assuming two independent motions; one at 
constant angle of attack, and the other at constant airspeed. Using. this assunp- 
tion, Eqs. (2+6) were separated and reduced as shown below: (negligible 


values from Table III are omitted) 





For Phugoid motion with < 2) 2958 =Q ; cos®~. = 1.0 
wo : Uo O 
S - Ay g u Ar 0 T 
Eqs (2-9) = 
-“% ~-S§ Q 0 Zn n 
Uo Uo 
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Then 


Eqs(2-10) | Xo g 








Eqs(2-11) -g Ly 
Uo 
u 
——— — 
7) CZ. KB = Zy 
u 7 Ww 





Oo 


For the Short Period Motion with u = 0 


Eqs(2-12) /S - Z, -S | “| o & J 
O 





- Uy (M,+SMy) S24-M,S g My, My] \7) 
Eqs (2-13) 0 - s 
ll a Mr g2 - M,6 
I S - Zw - § 


- Up (My +SMy) S2 - M.S 


ae Mr 
T S2 - (2, + Mg + UpMy)S+2Z,M, - U My 
Eqs(2-14) 
Pigamioc(2 
ce i ee 
ee eae ee 
7) S (Zw yf Mg ° UgMy,) one ZwMg UoM,, 


In this way the aircraft equations of motion reduce to second order 
Laplace equations which are extremely simple to handle. The time response 
to any type of input of either elevator or thrust can now be approximated by 


standard Laplace techniques. The response of airspeed to either thrust or 
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elevator will be controlled by the lightly damped (at this airspeed) phugoid 
characteristic equation while the angle of attack response to either forc- 
ing function will have the characteristic of the short period mode. 

At this point a further simplification can be made by examining the 
short period motion. By substituting the aerodynamic values in the above 
equations it can be seen that the short period mode has a damping ratio of 
. 35 with a time to one tenth of the initial disturbance.walue ofG..6iseconds 
indooeniod efdg yseconds. 

Thus it can be seen that the response to this disturbance: will. ap- 
proachcamiexponentially damped time response rather than an oscillatory 
type. This leads directly to the well known assumption that the aircraft 
moment of inertia in the pitch plane can be neglected. In the Laplace 
equations above it can be reflected by eliminating the S2 term from the 
matrix. : 


The equations then reduce to: 








Eqs (2-15) -. My) - qa 
“7 OS (= UGMy, - Mg) + ZyMg - USM, 
—_ a 
T = Vani eM, ier ec USN 
Which are of the form 
Eqs(2-16) eS ee 
S +A nl S +A 


lthis is more evident from the NACA nondimensional form of the equa- 


tions. The coefficient of the $2 term is then much smaller than the other 
terms. 


where A is 


Por a step input of elevator or thrust, Eqs(2-16) transform to the time do- 


main as 
Eqs(2-17) CC (t) al Cy aut _X (t) E Co ~At 
To aXe & Tt) “a \t - € 


Since A is usually large. (3.5 in this case) the exponentionaj term can be 
neglected. This indicates that the angle of attack response to either thrust 


or elevator will be approximately constant. In aerodynamic terms, 


Eqs(2-18) H(t) " My ~ Zn Mo/ U5 XK (t) = Mr 
T(t) ~ ZyM, - UpMy The) eZ, yhtge - igh, 


Substituting values gives, 


Eqs (2-19) A(t) _ _ A(t) _ _ -6 
my 7 7 )74 Ty) = 71:2 * 10° “rad/lb 


The value obtained for «(t)/ T) (t) is within 10% of the value deter- 
mined previously by analogue results. [9] The low value for “(t)/T(t) shows 
that the angle of attack response to thrust is negligible. The Power Com- 
pensator then can be used oT to control the airspeed perturbations .and 
elevator must be used to control the angle of attack perturbations. This is 
not to say that the airspeed variations will not affect the angle of attack. 
Any airspeed perturbations will cause the aircraft to seek a new equilibrium 
angle of attack. The angle of attack loop is necessary to offset the large 


rise in drag for positive excursions of angle of attack. 
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The assumption of two independent motions is now carried into the 
control system of Figure p by considering each loop separately. The equa- 
tions are further simplified by assuming all delay terms are equal to unity. 
The two simplified loops with their closed loop transfer functions are shown 


in Figure 4. These functions with aerodynamic terms substituted are shown 





below. 
Eqs (2-20) u _ - XK, (S + ky) 

ug SZ + KR - XS - GOA + RK, 
Eqs(2-21) _&% T Mr 


Ke S2 + (2, + M,)S + Z2yM, - UpMy - MKe 


Roots Locus plots for these functions can be easily estimated since 
they are only second order. Root Locus plots for these functions are shown 
in Figures 5 and 6 for the aerodynamic values listed in Table II, 

The gains can then be chosen from these Root Locus Plots. As a first 
approximation the gains can be chosen such that the phugoid roots are on 
the imaginary axis in close proximity to the short period roots. This will 
give a awenct of moderate overshoot with the damping of the short period 
roots. Note that the short period roots are changed only slightly from their 
basic aircraft position and that the gain of the angle of attack loop affects 
these roots only to a small degree. 

The angle of attack gain can be approximated by considering its de- 
sired function of equalizing the drag rise due to the angle of attack. This 


gives 


—— 


¥.2 = eat where Tc = Ky Ke 








Therefore Ke = “3 pv4sc be This gives an approximation of 24, 000 lb/ 
rad. The values of the other gains chosen were 
Ky = 300 - 500 lb/ft/sec 
ky = .05/sec 
The effect of the gains on the overall motion can be analized from 
the simplified Root Locus. Ky will mainly affect the time constant of the 
response. . Thus incredsing.K 2will increase the ‘speed of the response with 
a resulting increase in overshoot. The effect of K« is small but will increase 
the overshoot as itis increased. The effect of k,, is difficult to analize 
without a detailed study of the Root Locus. In general, increasing ky will 
increase the overshoot and can cause a second oscillatory mode to appear 
at about the same frequency as the short period mode. This would give an 
undesirable= response. 
The gains chosen as first approximations are ¢ 


K,, = 400 lb/ft/sec 


ky .05/sec 


Ke 


20,000 lb/rad 


Kuky = 2.0 lb/ft 


2.8 Co-ordinated Elevator Input 
In order to stabalize the aircraft motion at a steady value of airspeed, 
elevator must be used to control the flight path. This input will simulate 


the action of either an autopilot feedback or the action of the pilot. The 
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gain for this elevator/airspeed loop is obtained from the following analysis; 


V20=(eetou)2 y= “2.  Senecsaw 
BP CE 
then 2U du = y. - 2 dC; 
ZPCL 
or 2uU — A CLaet 
ZPCr 


therefore 


ee ae 
puc;? 


For the given flight condition 
u(kts) =-140.0“<(rad) = - 6.32(deg) 
From Eqs 19, oc = - 1,747) 


Then u(kts) = 11.07) (deg) 


3 Analogue Simulation. 

In order to check the validity of the foregoing analysis, the full 
equations (Eqs. 2-4) were solved on an analogue computer and the gains 
of the Power Compensator "optimized" for horizontal and vertical gust in- 
puts. A Donner 3100 computer was used in the analysis. The analogue 
wiring diagrams are shown in Figures 7 and 8. All analogue program infor- 
mation is given in the Appendix. The values given for the flight condition 
in Table III were used in the program. Figures 9, 10, and 11 show the 


response of the basic airframe to :horizontal and vertical gusts and to a 
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step elevator input. The motion shows the phugoid response with a period 
of 28 seconds which compares with the calculated value of 30 seconds. The 
short period motion is difficult to see because of the heavy damping but the 
period is less than 10 seconds. Figure 11 clearly shows the relation be- 
tween elevator input and angle of attack. The angle of attack rapidly takes 
on a small phugoid oscillation about a constant value of - 1.8 deg for an 


elevator input of 1 deg. This compares well with Eq(2-19). 


3.1 Basic Airframe 
Figures 12, 13, and 14 show the response of the basic airframe to 
sinusoidal horizontal gusts with periods close to those of the short period 
and phugoid motion. Figures 12 and 13 show that the response to the short 
period gust is of small magnitude and composed of the short period motion 
superimposed on the longer period phugoid. As the period of the phugoid 
mode is reached, Figure 14 shows that the response becomes divergent. The 


same is true for vertical gust inputs shown in Figures 15 and 16. 


3.2 Basic Airframe with Power Compensator 
In order to check the validity of the gain analysis given in 2.7, the 
response of the airspeed motion for various values of loop gains is given in 
Figures 17 through 20. Figure 17 shows that K,, affects the damping of the 
motion and is the main controlling gain. Figure 18 shows the decrease of 
damping with increasing integral gain, Kuki Figures 19 and 20 show the 


small effect of the angle of attack gain, Kx, on both horizontal and vertical 
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gusts. Since the power compensator will have the most effect on the air- 
speed loop, the airspeed response was used to choose the "best" gains. 
Figures 21 through 25 show the response of all the important para- 
meters to various gains for the most unfavorable gust condition; - 5 kt 
horizontal and + 5 kt vertical gust. Important aspects of these Figures are 


summarized below. 





TABLE V 
P _. ANALOGUE DATA 


Bic Overshoot within 1 kt . required 10 sec 





21 2000 lbs 
22 : 1800 
23 - 2800 
24 : 3500 
29 : 1400 


The gains from Figure 21 were considered as the most favorable. The over- 
shoot is somewhat misleading since it will be reduced considerably by the 
pilots application of elevator to correct the flight errors. However, this 
overshoot could be reduced by series compensation. 

The gains of Figure 21 were used for all the remaining runs containing 


the power compensator and are listed below. 


K,, = 400 1b/ft/sec 
Kuku = 2.0 lb/ft 
K« = 10,000 Ilb/rad 


Figures 26 through 30 show the response of the airframe and power compen- 


sator to sinusoidal gust inputs. They show that the compensator eliminates 
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the phugoid divérgence but Figure 29 indicates that the system tends to 
magnify the — for vertical gusts with a period of 9.6 seconds. How- 
* ever this period should be long enough to be controlled by the pilot and the 
airspeed perturbation is not serious. Figures 31 and 32 show the response 
to an airspeed command with coordinated elevator input as discussed in 2.8. 
These Figures show that the overshoot is reduced to 26% and the glide slope 
variation to within .5 degrees. 

4, Root Locus Analysis. 

In order to compare the simplified Root Locus plots with those for the 
complete system, the characteristic equation for the complete system was 
evaluated from Eq(2-8). This is shown in the Appendix. Root Locus plots 
of this polynomial were then developed for several values of loop gains. 
These plots are shown in Figures 33 and 34. 

Se Conclusions. 

1. The airspeed/angle of attack power compensator would be a feas- 
ible system to employ especially where the use of a full auto-pilot is anti- 
cipated. 

2. The simplified method of setting loop gains by considering the 
control loops as separate functions of the phugoid and short period motions 
will give valid results. 

3. The response of the system can be approximated by considering 
the separate root locus plots of each loop. 


4. The system may have to be improved because of the excessive 


Zo 


overshoot response of airspeed. This could be done by series compensation. 

S. The effect of the gains on the Root Locus of the complete system 
compares qualitatively to the effect of the gains on the two loops consider- 
ed separately. 

6. The power compensator will have negligible effect on the short 
period motion and can be thought of as an airspeed (phugoid) damping device. 
This will be true whenever the thrust axis is close to the body axis giving 
small values for My and Lips : 

Op Recommendations for Future Study. 

1. The effects of the non-linear aspects of oS. on the pitch motion 
feedback gain requirements for all three power compensators studied thus 
far-should be the next field of study. 

2. An investigation to determine which of the three systems studied 
thus far would give the best results when combined with an elevator auto- 
pilot:would indicate an optimum combination for a complete system. 

3. A complete analysis of a flight control system should include the 
effects of turbulentair gusts and aeroelasticity. To do this type of an analy- 
sis by analogue simulation would require an excessive number of compon- 
ents. Since the problem of response to random frequencies is best analyzed 
by frequency response methods, a digital computer program would give the 
most convenient approach. This type of a problem could also be analyzed 
ona "hybrid computer system" by using the analogue to solve the equations 


and the digital to introduce the non-linearities. 
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FIGURE 3 
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FIGURE 5 
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7s APPENDIX 
7.1 Digital Computer Programs 

The digital computer programs are designed to compute the character- 
istic equation and the transfer functions for elevator and thrust inputs. By 
utilizing the proper combination of subprograms either wind tunnel data and 
airframe dimensions or reported nondimensional derivatives may be used in 
these programs. The subprograms developed here include only the longitud- 
inal motion; however, any number of other subprograms may be added with- 
out revision. The program language is defined in Table WL. All input data 
must be in the foot-pound-radian system unless otherwise noted. Standard 
programing proceduress for FORTRAN 63 are used. [1] The subprograms and 
a sample print out are shown in Figure 35. 

Several of the subprograms require the use of another subprogram to 
compute the roots of the generated polynomials. These subroutines are of 
a standard type but the usage will vary depending on the installation. The 
subroutine used in these programs to compute the roots is RTPLSUB which 
is a library routine available at the Computer Facility of the U. S. Naval 
Postgraduate School. However, any subroutine designed to compute-roots 


of polynomials could be used by rewritting the call statements. 


7.1.1 Main Program - AEROFUN 
The main program may be of any desired format so as to utilize the 
necessary subprograms and print the results. The main program used here 


is dimensioned to take up to 50 sets of data for different flight conditions. 


fie 


The variable NN is used to indicate the number of sets of data for analysis. 


7.1.2 Subprogram COMPTDER 
This subprogram computes those aerodynamic derivatives which must 
normally be approximated from the wind tunnel data. The equations utilized 
for this program are listed in Table VII below. [5] 
TABLE VII 


EQUATIONS UTILIZED IN COMPUTER PROGRAM COMPTDER 





ye 27 CoS “A Cy 
ig ese se 
bc * [ere ney IRS oP 
Si sa 
Chow = MMP 2 (Xa — Kal Gy Cu =eeacuuee Cae 


Crm gw a Chaya 5] 























ve Arecosh | si AR+acosf \| - 16 
Bc 8 AR — | Rt6cosf | ee 
F $e\" 
_ t = woke 
oe ee Cingg = ~-O4 aire 
Cig = Cilee + Or. Cwgq im Craw + Cmae = Cra. 
S. 
Cu se | a — 
€ 
Cm = -22t a 7 & 
CaceZeez Maiice: TelaC Cmge = ELA + BC 
Mm ¢ge oo ce ft bey wigs Ge Sw Lal 
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INPUTS: A A Cr, x: Cla, r) AR, C, Swe hy, ih), ’ St. Fy A, B, Xacs Xcgi 
O€ 
Li, OK ® 


NOTE: If CL, is based on the wing reference area, assign Sy and 5} values 
of one. 


o 


Se mY either be an input or it may be computed in the program. If 
it is to be an input, assign the value to DEDAI. If it is to be com- 
puted, assign zero to DEDAI, 


Oe€ 
OUTPUT: Crq: CLge + Cmar Cmges Cmser Smee SE. 


7.1.3 Subprogram ENGDIMEN 
This subprogram computes the English dimensional derivatives from 
the NACA non dimensional derivatives. The relationships used in this sub- 
program are given in Table II. 
INPUTS: U.P. W, Sy. Be G prs 2qe lye Cyr Cher Cry Cig: “aie : 
Cp, Spx: Cpu Chq CDice Some: “me Cras Cmee OT .- 
COTS: Ka Xyg, Xu Xy » Xp, Zy, Swe Zqe Sq» Sp, Mye My, My, Ma. 


My. Mrz, g cos8,, g sin 85: 


7.1.4 Subprogram STAB 
This subroutine computes the coefficients of the characteristic equa- 
tion, the roots of the equation and the characteristics of the short period 
and long period motions. The equations used in this program are given in 
7.2. 


OPUTS: U,. Gg. Xye Kaye Zaye Zaye SQe he Migs May My. 
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OUTPUTS: (as defined in Table VI) U, V, CONV, WNLP, WNSP, ZLP,, 
Z5P/ TEP, ieee. Hor, TTLP, tisk. 
NOTES: Requires a subprogram to evaluate the roots of the generated poly- 
nominal. The coefficients of polynominal are stored in the dimen- 
. o.. » gloned variable C. The order of the polynomial is assigned to the 
variable N. The outputs of the called subroutine (the roots of the 
polynominal) are assigned to U and V and are therefore preserved 


for subsequent calculations. 


7.1.5 Subprogram FORCEFUN 

This subprogram computes and prints the coefficients and the roots 
of the polynominal numerators of the transfer functions HST): u/%) ,X/T, 
u/T. The equations utilized in the program are given in 7.2. 
INPUTS: Uo, 89, Xu, Xw. XT. Zur ws 2qs 2n, My, My, My, Mg, My, 

My . 

OUTPUTS: None. 
NOTES: Requires a subprogram to evaluate the roots of the generated poly- 


nomial. The operation is the same as for STAB except the output 


is printed internally. 


7.2 Derivation of Airframe Transfer Functions 
The airframe transfer functions are derived from the equations of 


motion written in Laplace Matrix form given previously as Eqs 2-6. The 
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TABLE VI 
DEFINITION OF SYMBOLS 


Program 


Mnenonic Parameter 
AR FR Aspect Ratio 
B B (Iw, slug-ft“) 
CF Cn c chord 
CD Cy 
CDA Cy. 
CDDE Cy. 
CDQ CD, 
CDU Cay 
CL Cy 
CLA Cr x 
CLAT Cr ets 
CLDE oe 
CLG Chg 
CLU Cru 
CMA Ck 
CMAD Ch. 
CMDD C..¢ 
CMDE C150 
CMQ Crng 
CT Cn 
CONV Convergence indicator for roots 
output of STAB only 
D , Zp Perpendicular distance between 


the cg and thrust axis. Pos- 
itive for thrust axis above cge 


1S 


DHDA, DEDAI 
ET 
GCTHETA 
GSTHETA 
ET 


ST 
sw,s 
TAU 
THETAD 
TLP 
TSP 
THLP 
THSP 
TTLP 
TTSP 


U,U0 


W 
WNLP 


WNSP 


TABLE VI(CONT) 


Ne 
aA 


YN, effectiveness of horizontal tail 

g cose 

g@ sine 

h, height of root chord position 
of horizontal tail above wing 
root chord. 


Number of sets of data =~ Main 
Program only 


P (slugs/ft) 





St 
S, or reference area 
"~ oat 
Vda 
®, (deg) 


Period of FPhugoid Motion 
Period of Short Period Motion 
Phugoid time to 4 

Short Period tim to 4 
Phugoid time to el 
Short Period time to el 


“A taper ratio 


Ug wveference or steady state 
airspeed 


Output of real parts of roots 
STAB only 


Output of imaginary parts of 
roots = STAB only 


weight 
Wy phugoid 


Wy, Short period 


7 


TABLE VI (CONT) 


AE TA | Ky 
XLAMBDA IN. sweep (deg) 
XL 1, distance from cg to aero~ 
dynamic center of horizontal 
tail 
UU M, 
XMETA My 
ALIQ Mg 
XT . Mm 
XUeW My 
WD Me. 
XU ne 
XQ Xx, 
XT Xp 
xi Xe 
ZE yw engle between thrust line and 
body axis, positive if upward 
(deg) 
ZET Zy 
ZQ Za 
7a Zm 
ZU Z 
zu 2 
Ze? ° of phugoid 
ZSP “ of short period 
A&B COMPTDER only - values from 
Perkins & Hage, Fig 10-10. 
SAC Xac position of aerodynamic center 
of wing in % MAC 
XCG Xeg position of center of gravity 
in % MAC 
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characteristic equation is the expansion of the main determinant as shown 








below. 
S-X, — Uae, g cos Oz 
= ie & és - |e | g SIN O, 
CH = U, a ae U. +i5 + UL 
2 
—M, — -U[Mu - SM) S°— MeS 


CH= S'T+HAS' + BS' + CS + D 





A = — Xu —-Z, SIME —M, Lo t U, 
5 = Xa Les gi Mg|Xut Tm) — a Zs a MIE. mn U,| a Mol Xul Le a U,] 
+ 9 SIN e,] 


C= Myl-Xy Zw + Xw Zu) + Mulg cos Go— Xw(Uo + Za 
+ Me(Zug cosQ,- Xu sin Og) + Mw] Xu(Zq + U,] +9 SIN e,] 


De g cos Oo(My Zu - 7a + g sin O,(M Xa — Xu Mw] 


The transfer function relating airspeed to a thrust input perturbatiom 











is then: 
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The other transfer functions follow in a like manner. 
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7.3 Evaluation of the Closed Loop Transfer Function. 


The closed loop transfer function was derived as Eq. (2-8) 


do ff Te au 
Uc ibe 2 2. Oy, See ES 
Uc Te Te T Ke Té Tf T 


ae U 
eee 
Uc i¢ t¢ vw _ 1c X 

ue T Ke T 


Substituting the transfer functions defined in Table IV with 


ee Sy Sell 
T CH T Gr 
gives 
KulS + ky) uN 
We, (a a eee 
Us 1 + Ru(s + ky) UN _ xX oN 
S CH “ CH 
therefore 
— K (S + k,,) (UN) 
Me S(CH) + K, (5 + ky) (UN) - K« (*N) 


The denominator of this expression is a 5th order polynomial shown below. 


So + c,s4 “ C,8° a3 C,82 t C45 + Cec where, using 


the same notation as 7.2 


Co = B+ KyXpAyy + KykuXp - Ke Aang 


C3 = © + KyXpByp + KykyXpAypr - Kae Ac qBe 7 
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Cg = D+ KyXpCyup + KykyXpByp - Ke Anat 


Cs = KykukeCur 
Then 
u_ _ _Xaky(S + ky) (63 + AyrS? + BytS + Cur) 


From this expression the root locus was obtained by digital computer means. 
The integral gain produces one variable zero. The other zeroes of the func- 
tion are invariant with changes of the gains. The root locus plots obtained 
are shown in Figures 33 and 34. The most desirablee gains can now be 
chosen from these plots using standard pole-zero methods. [2} However, 
this type of an analysis would have to be made for each loop involved to .. 
arrive at the best combination of gains for both the airspeed and the angle 


of attack response to gusts. 


7.4 Analogue Simulation 


Substituting the aerodynamic values into Eqs 2-4 gives 


u = -.060U°- 3.32% - 31.850 +..00145T - 1.647) 
K& = -,00l13u - .4265< +6-.019350 - .08157) 
6 = 1.85xl07~4u - 1.14% - .0415& - .3386 - 4.55x1l0-6T 


= "2 20u) 
These equations were then scaled using the relations defined in Table 


VIII. Real time was used. Using this scalirnigs the voltage form of the equa- 


tions becomes 


G = -.10lu - .O116< - .1110 + .292T - .0067) 
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K= -.544u - .4265% + 6 - .019358 - .08157) 
@ = .03llu - .398% - .0145%.- .12e06 - .7857) 
- .091T 
The power compensator equation was scaled as (Eq 2-5 without 


delays) 





Seal ee dl al 118.9 “é 
The delays were simulated by the general relationship shown below 


were 9, is the delayed signal and 9, is the input signal. 


8 _ ] = . 
If 8; = Te jac: then Bo = 1/T Qj 1/T 





or 8, = pre = 70 oe 


which gives the simulation circuit 





The values of 1/7 can be obtained by either potentiomenter settings - by 

adjusting the input gains. a 
The schematic diagram for the wiring circuits is given in Figures 7 

and 8. The number on each input is the gain of the.input where the gain is 


as defined below: 


For a summing amplifier For an integrator 


R¢ 
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R > 
—_ — — Slee 
Gain Ri Gain Ric 


The potentiometer settings are given in Table IX, In order to keep 
the witings diagram as simple as possible, feedback circuits are indicated 
by the use of numbered potentiometers or numbered units shown dashed. 
The squares indicate recorder channels. SW 1 and SW 2 are internal double 
throw-double pole switches. NC refers to the "normally closed" position 
of SW 2. Step gusts were simulated by means of initial conditions on the 
airspeed and angle of attack integrators. Sinusoidal gust inputs were 


obtained by means of a frequency oscillator. 
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TABLE: VIII 
SCALED VARIABLES 


S| 
I 


50 ft/secé 
8.3 ft/sec 
50 


5 ft/sec/volt, 


0843 ft/sec/volt 
~5 kts/volt 


0001715 rad/volt 
ol deg/volt 


cls 
uit 


kts 


0175 rad 
LO deg 


0005 rad/sece/volt 


20017L5 rad/volt 
ol deg/volt 


5 rad/sece 


ol7L5 rad 
LO des 


HI] Diny DI] HA 


1000 ibs 100 Lb/volt 





98 


Pot # 
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Gee ie a 


TABLE IX 


POTENTIONETER SETTINGS 


Sevting 


0292 
0006 
OL 
00116 
der 
0101/.0116 
20815 
5a72 
wee 
0015L 
027242133 
0091 
0393 
0031 
099 
eOUs 
mils 
2031/0 
Various 
0/ +199 
hy 

K 


oy 


ieioyb gram 


Parameter 


= 


o ce 5 


Up/ Ke (Gusts) 


Up/K,, (Gusts) 
© 

Cc 

Up/K, (Gusts) 


Us 


fe 


Gain 


PrP FP PP 


im 


tt 


2/2 


mM FY FY FSF Ff Nf FP YS 


ee oe oe 
oO 


31&32 


‘ 
33&3h. 


ee 
36 
38 


TABLE IX (CONT) 


5 

olt35 
Various 
Various 


=) 


T, (Throttle Servo) 
Te (Engine Delay) 
n 

Vc 

VY delay 
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